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Abstract 

This paper addresses a specific issue of time accuracy in the calculations of external 
aerodynamic problems. The class of problems that is discussed here consists of inviscid 
compressible subsonic flows. These problems are inherently governed by a convective equa- 
tion. This is readily seen by linearizing the Euler equation which results in a convective 
wave equation for the pressure. A key mathematical issue that is not well understood in 
literature for these problems is the long time behavior of the solution. This is an important 
aspect if one desires transient calculations of problems governed by the Euler equations or 
its derivatives such as the small disturbance equations or the potential formulations of the 
gust problem. In particular, difficulties arise for two-dimensional problems. In two dimen- 
sions the time decay rate of solutions of the wave equation is known to be rather slow. This 
applies to the convective wave equation as well. The consequences are rather severe if one 
focuses on the time accuracy of solutions of problems governed by the Euler equations. 

In concert with the above mentioned problem, exterior flows require proper modeling 
of boundary conditions. In particular, computations of these flow problems require trun- 
cation of infinite regions into finite regions with the aid of artificial boundaries. On these 
boundaries one must impose boundary conditions that are consistent with the physics as 
well as guarantee consistency with the original problem posed in the unbounded region. 
Moreover, these boundary conditions must have accuracy in time as well as space. Some of 
the well-known procedures do address the issues of spatial accuracy and have remedy for 
these conditions. Unfortunately, these procedures do not address the time accuracies, which 
are crucial for the transient problems. Our treatment is discussed in detail and examples 
are presented to verify the results. 


'Research for the author was supported by the National Science Foundation, DMS-8921 189 and by ICOMP, 
NASA Lewis Research Center (work funded under NASA Cooperative Agreement NCC3-233). 
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1 Introduction 


This paper deals with unanswered aspects of long time behavior of unsteady flow computations 
in two dimensions. Many of the available codes such as the FLO series deal with steady flows 
and their accuracies are well tested. State of the art grid generations and even complicated 
geometrical structures are included in these codes. While these developments reflect the ma- 
turity of computational fluid dynamics, the goal of the field is to understand aerodynamics in 
its entirety which includes the unsteady phenomena. It is important to understand this phe- 
nomena if one is to achieve successful calculations on biproducts such as far field acoustics. A 
common tendency is to use steady fluid dynamics codes for unsteady calculations. Such cal- 
culations yield either unsatisfactory resolution of time accuracies or break down in finite time. 
The classical problem of sub sonic- transonic compressible flows past a cylinder is an unresolved 
issue for this very reason. For example one can see a documentation of such computations in 
[15], The key reason why this is difficult is due to lack of understanding in the modeling of 
far field boundary conditions that naturally plays an important if not an essential role in these 
computations. 

It is an open issue whether one can predict solutions of unsteady Euler or Navier-Stokes 
equations in a satisfactory manner. A recent paper by [2] deeds with numerical computations of 
viscous flows past a circular cylinder. The modeling of required far field approximate pressure 
condition was obtained using a theory proposed by [1] which is accurate for high frequencies. 
This limits the solutions to be accurate for short times. They also proposed approximate 
boundary conditions using characteristic variables which works better than the above said 
pressure boundary condition. It is not clear if these conditions mantain time accuracy of the 
solution and it needs to be addressed. This concern forms the basis of this paper. It should be 
emphasized the workers realized similar difficulties ahead of time and concentrated on linearized 
Euler equations and in the frequency domain. Namely, the work of Goldstein [12] stands out in 
the literature. Corrections to this model were proposed by Atassi [13] and subsequently he and 
his coworkers (e.g. see [14]) succeeded in predicting unsteady aerodynamic calculations through 
numerical procedures. These formulations and calculations are done through potential methods 
and the resulting equations are nonconstant coefficient second order linear partial differential 
equations. Calculations are accurate within the linearized limits. In certain cases, such as flat 
plate air foils, the theory behind these equations is well understood. As a result, these potential 


2 


formulations have firm foundations. Recently several authors are beginning to solve linearized 
Euler equations directly. It is not clear if there are any real advantages to these procedures. 
It appears that the mathematical foundations do need careful work. Moreover, associated far 
field boundary conditions consistent with the linearized equations require a careful study. This 
is the other aspect this paper addresses. 

Two-dimensional problems governed by any wave equation, unlike their one- or three- 
dimensional counterparts, suffer severe slow decay rates under certain conditions. Whether 
the problem is posed in the frequency domain or in the time domain, the problem of slow decay 
rates occurs. For example, for a very low frequency problem it would appear from a concep- 
tual point of view that the terms multiplied by the frequencies could be dropped to consider a 
problem in the absence of those corresponding terms. A typical situation is the problem gov- 
erned by the Helmholtz equation, where reduction is the Laplace equation. Similar arguments 
hold for the time dependent wave equation. If in the long time there is a steady state, the 
analogous situation is to consider the steady state part. These discussions imply that there is 
interplay between low frequencies and long time behavior. The researchers who looked at these 
questions considered only isolated issues. In this paper an attempt has been made to connect 
the theoretical issues behind these problems to aerodynamic problems. 

In section two we discuss the behavior of the solution of reduced wave equations in exterior 
regions. In this section we point out the potential problems in considering the zero frequency 
limit or the steady state problems. Moreover, we show the proper form of the boundary 
operators one should consider in the prescription of radiation boundary conditions. In section 
three we consider a true time dependent model in which the corresponding corrections to the 
radiation boundary operator are discussed. Here we include available well-known procedures 
that are commonly used in the computational fluid dynamic literature to show these corrections 
in perspective. Numerical solutions are shown to establish the validity of the theory. 

Finally in section four we provide extension of the theory proposed in section 3 for the treat- 
ment in the far field. This treatment should hold for the Euler equations and their linearized 
versions. The corresponding numerical treatment will be discussed elsewhere. 

2 Low frequency behavior of the reduced wave equation 

Here we show a potential problem that arises with the reduced wave equation. In particular 
we shall discuss the low frequency behavior that arises from problems with Dirichlet boundary 
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conditions. The reason for discussing such behavior is to establish the slow decay rates of the 
two dimensional wave equation. These decay rates play an important role in the construction 
of approximate boundary operators that are vital to the constructions of solutions. At the 
outset we want to emphasize that these decay rates and associated problems only apply to the 
two-dimensional wave and wave-like equations. 

We shall begin with the two-dimensional problem that will motivate the decay rates in 
question. We consider a classical soft obstacle problem in two dimensions. In particular we 
consider the following boundary value problem (P^): 

A u + k 2 u — 0 in fl e C 7Z 2 

u - g k on T (1) 

lim (u T — iku) = 0 r = \fx 2 -\-y 2 

r— kx> V 

where Q e is the exterior to a bounded region f l and T is its boundary. For the soft obstacle 
problem where tw is the incident field. 

The key question we address is what happens to the problem as k — > 0? To discuss this 
issue we consider an appropriate zero frequency problem. This problem is constructed in such a 
way that it preserves the zero frequency behavior on the equation and the boundary condition 
on T. However, the radiation condition at infinity is replaced by a boundedness condition at 
infinity and is as follows: 

Av = 0 in fl c C TZ 2 

v = g 0 on T (2) 

\v\ < M or v — ► A log |x| as |xj = r -» oo. 

We call this problem to be (Po). It has been shown in the literature [10], [8] that in the 
limit as k — ► 0 the solution of (P^) approaches the solution of (Po). However, this limit process 
is a slow one and we want to make a case for it as it has everything to do with constructions 
of boundary conditions for numerical simulations. Before we introduce the results let us focus 
our attention on the structure of v mc = —go. The choice of Vi nc depends on Ui nc . Examples are 
as follows: 

a) ^inc = e%kx (plane wave) 

Then the corresponding zero frequency limit is readily seen and is t>j nc = 1 


4 


b) Uinc = -|^o 1} ( fc l* “ *o|) (singular source) 

Then the corresponding zero frequency limit is taken as the principal leading spatial singularity 
Vine = log | x — Xq\. A discussion of this choice is seen in [7] and is motivated by several issues 

including an obvious one: the asymptotic behavior of the Hankel function for small arguments. 
Now we state some of the key results: 
a) (Crude) : u — + v as k — ► 0 


b) (Sharp) : u = v + 0(^) 


c) (Sharper): u = v + + 0(k 2 (lnk) 2 ) 

The results in a) are due to an old work of MacCamy [8]. The results in b) are later 
observed by many others, namely Muravei [10], Werner [16] and Kress [17] which indicate the 
slow decay rates. The result c) is due to the author and his coworkers and the results appear in 
[9], [6] and [7]. To motivate the implication of these results first let us demonstrate the effect 
of these decays numerically. For this purpose we have taken a plane wave incident with unit 
amplitude. It is easily seen from problem (Po) that the solution is v s = — 1 and the total field 
v — Vi nc + v a = 0. This solution is compared with the total solution for the nonzero frequency 
problem (P*.) with u = Ui nc + u s . Figure 1 shows magnitude of u for the frequencies 10~ 2 , 
10“ 3 and 10~ 4 respectively. Clearly one can see that even for these very low frequencies the 
solution does not approach the correct low frequency limit of zero rapidly. The frequency or 
the wave number approaching zero corresponds to the long time limit of the solution to the 
corresponding time dependent wave propagation problem. In particular for moderate times or 
frequency one cannot expect to use any properties of the steady state problem and in this case 
the solution of the zero frequency problem (Po)» This message has serious impacts on unsteady 
flow computations. 

In light of the above discussion, we would like to outline a correction to this zero frequency 
computation. These ideas will later be transferred to an appropriate time domain problem and 
at the end we shall propose ideas as to how to extend the theory to the Euler equations. In 
particular we outline a theory that appears in [6], [7] to show the sharp result b) and the sharper 
result c). The construction of these results involve potential theory. We use potential theoretic 
representation for the solution. Specifically we use single layer potential representations. Let 
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(3) 


u(x) be the total field in ( P fc) , i.e., 

u(z) = U inc (x) + U,(x). 

Then, for the problem (Pk) the representation is as follows: 

u(x) = J^cr(y)Gk(x,y) dsg +u inc (x), x e (4) 

where 

G k (x,y) = ~H { 0 1 \k\x-x o|) (5) 

is the free space Green’s function. The scattered field satisfies the Sommer feld’s radiation 
condition due to this Green’s function. One obtains a boundary integral equation for the 
determination of <r(y) that is defined on the boundary F using the boundary condition on T. 

Jr<r(y)G k (x,y) ds# + u inc (x) = 0 ,x G T (6) 

A similar construction with the free space Green’s function for the Laplace’s equation is made 
for the problem Pq noting that any representation of the solution is subject to an arbitrary 
constant. 

l>(;c) — J Gq(z!, 2/) -f- Vtnc(®) “he, x E (7) 

with 

j^{y)ds g =- 1 (8) 

where 

G 0 (*,y) = ^(l°g|* -»|). (9) 

The condition described by equation (8) takes care of the arbitrariness of the constant c and 
the boundedness of the solution at infinity. Again demanding the boundary condition on T we 
arrive at the boundary integral equation 

J^n{y)G 0 (x,y) dsy + v inc (x) + c = 0, x e T. (10) 

Thus the complete solution of (P 0 ) is determined through solving the equations of the pair of 
equations (8) and (10) simultaneously. The key idea of the low frequency correction relies on 
relating cr and fi with appropriate low frequency corrections. Several of the quantities that arise 


(ii) 


can be related to the asymptotic formula for small values of z 

~ M = Yr (Mf) +7) - | + 0(1 J lnz) 
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where 


7 = .5772156649 • (12) 

is the Euler constant. Using this expansion in equation (6) we obtain a related integral equation 
in the following form 

J <r (j7) G 0 (x, y) ds# + j3(J^ a (y) dsy + 1) + Vinc(z) = 0 ( 13 ) 

where 

13 = h (log( f > + 7) ~ *■ (14) 

Successive approximation by considering equation (4) and (13) yields ([6]) 

a —a + 0(k 2 Ink). (15) 

Further it has been shown in [6] that 

<r (x) - n(x) + C J ^1 , x e r (16) 

P - Co 

where / 0 ( x) and the real constant c 0 satisfy the following equations : 

j fo{y) Go{x, y) dsy + Co = 0 (17) 

Jfo(y)dSg = 1. (18) 

Use of this result leads to the results in c). In particular two low frequency corrections can 
be constructed. One is for k\x\ small (the results in [6]) and another one that is proposed in 
the recent paper [7] for k small globally. The emphasis of the latter corrections are to include 
the low frequency corrections both in the near field as well as the far field. Those corrections 
are as follows 

u> (x) = J O' (y)Go(x,y)dsg + (y)dsy + 1) + tw(iT) (19) 

which is the one that appears in [6]. 

ty(f) = J a {y) (?*(£, y) dsg + iw(£) (20) 

This is the improved result that is valid both in the near field and in the far field. To illustrate 
the effectiveness of these approximations we solved for u, the full solution of Helmholtz equation, 
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and the approximations w and w for the case of a circular obstacle and a plane wave. All the 
results were calculated numerically. Figures 2 and 3 show the effect of these corrections as 
predicted by the theory for reduced frequencies k = 10“ 2 and k ~ 10 -4 respectively. 

The important information that is related to the computations of exterior aerodynamic 
calculations is as follows: Suppose one wants to compute solutions within a truncated region 
with the aid of an artificial boundary. Then for the well-posedness of the problem within the 
truncated region one must have normal derivatives prescribed on the boundary in an appro- 
priate form. For example one can realize this form by considering the Sommerfeld’s radiation 
condition. The operator that we propose is 

>,+ 0(*W) 2 ) (21) 

A similar relation holds if one is interested in specifying the normal derivative of the scattered 
field only. In this case u will be replaced by u 9 and v will be replaced by v s . The importance of 
this condition is that in the low frequency limits it is important to correct the zero frequency 
operator with the logarithmic terms which play a significant role. This will be evident in the 
next section where we shall deal with a true time dependent situation and the structure of 
these logarithmic corrections in the frequency translates to logarithmic corrections in time in 
the artificial boundary operators. These constructions were made from a computational point 
of view but are consistent with the abstract theoretical results posed in [10]. 

A point of comfort and warning is that these corrections are not necessary in the case of 
a hard obstacle. That is, the Dirichlet condition on T is replaced with a Neumann condition 
= ~ c ' ft can be shown that the decay rates are much faster for this condition and 
does not involve any logarithmic corrections. Potential formulations that were introduced by 
Goldstein [12] and Atassi [13] fall into this category. As a result the computational results 
presented in [14] do behave well in the low frequency limits. However, issues dealing with the 
time dependent problems, in particular problems governed by the Euler equations, are in need 
of further theoretical backing. We show a preliminary study in the next section. 

3 Implications to a True Time Dependent Model 

In this section we motivate a time dependent model problem that is used to indicate potential 
troubles in two dimensions. Specifically, we shall focus our attention on construction of ab- 
sorbing boundary conditions on artificial boundaries. The model problem that we propose is 
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governed by the wave equation exterior to a circle of radius 7*0 with a time dependent data on ro 
which asymptotes to a constant value as t -► 00 . For computational purposes we truncate the 
infinite region by an artificial boundary at v — R. At this boundary there is an exact condition 
(in the mathematical sense) that takes care of the effects of extending the region to infinity. 


The model problem is as follows: 

u tt - u rr + -u r , r 0 < r < R (22) 

r 

u(r 0 ,t) = 1 + g(t) (23) 

^ = !F(u) at r - R. (24) 

ti(r, 0) = u t (r, 0) = 0 tq < r < R (25) 

g(t) = 0{^) as t -* 00 (26) 

This model is designed to have the long time behavior of a solution u — * 0 as t — » oo. The 


operator T(u) is the exact operator known only theoretically. It has a well defined structure 
in the transformed domain which we discuss now. The initial values suggest that it is easier to 
work with the Laplace transforms. Doing so one obtains 


s 2 u = u rr + -Ur, r Q < r < R (27) 

r 

u(r;s) = g (28) 

») = ?(B) (29) 

The operator J- in the transformed space is readily verified by considering the solution of 
equation (2.27) in the exterior region and is as follows: 


JF(u) = 


sKq(sR)_ 
Kq(sR) U 


(30) 


Available well-known boundary condition procedures rely on obtaining approximations of the 
symbol or the ratios of the modified Bessel fucnctions. Here we provide some of these approxi- 
mations and operators. Let B be the approximation of J 7 . Then 
1. Engquist and Halpern [11] 


B = -s 


(31) 


9 


which is the crudest approximation of T. Then the approximate boundary operator in the 
transformed space is 


du — 
* * 


Inversion of the transform results in a time domain operator 

0 


du du 
dr + dt 


(32) 


(33) 


The important part to consider here is the error in the approximation (2.31) using asymptotic 
properties of the modified Bessel functions for small and large arguments. This in turn will 
suggest the order of the error in real time. Doing so we obtain 


7- B = 


i 

27? 


as s 


as s 


00 

0 


(34) 


2. Bayliss and Turkel [4] first order 

Analogous list in this are boundary operators in the transform domain 


B = -s - 


2 R 


the boundary operator in the time domain 


du du 1 
dr + ~dt + 2R U = ° 


and the errors behave like 


T-B = < 


I" O(^) as s — ► oo 

“275 as s -» 0 


3. Bayliss and Turkel [4] second order 

Here the boundary operator in the transform domain is 

c 2 I 3i | 3 

B = 5 Z 2R z sm 


5 + a 


and its translation in the time domain is 


. d d 5 w d d 1 . 

( d~r + ^t + W )( ¥, + at + 2R ) ' i - t, 


The errors in this case are 


f O(^) as s -> oo 


T - B = < 


m as s -» 0 


(35) 

(36) 

(37) 

(38) 

(39) 

(40) 
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4. Hagstrom [5] (part of current study) 

The boundary operator in the transform domain is 


B = - 


S Sj + 2R + Q 

s + 8 


and its translation in the time domain is 


( 41 ) 


. d c . , d 1 c .du 

{ Jt + {)( + ( a< + 2R + S) at ~ 0 


Finally the errors behave like 


T -B 


= 1 °(i) 


O(th) 


as s 
as s 


00 

0 


(42) 


(43) 


We propose to correct the errors in the low frequency limit, i.e., the long time behavior in 
constructions of these operators. This is the punch line of this paper and again we emphasize 
that this idea is crucial if one were to compute directly time dependent problems. For this 
purpose let us observe that if one were to include the next order correction to the asymptotics 
for small s, we obtain the following properties: 


( “(^+ 2 ^)^ as s — ► oo 


U r ~ < 


as S — ► 0 


(44) 


where j is the Euler’s constant. Now consider the case where 


u — ► 1 

which has the behavior in the transform domain 

1 

u — ► ” 


as t —+ oo 


as s — * 0 


In particular 


1 


tX r ~ — 


as s -+ 0. 


Rs In (Re7%) 

To obtain the corrected behavior in the time domain we look for a function G(t) so that 

1 


C{G{t)u) = - 


Rs ki(Rei |) 


as s — * 0. 


(45) 

(46) 

(47) 

(48) 
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Construction of such a function will be found in [3]. The key result here is Proposition: 


G(t) = 


( 49 ) 


i2[ln(^)+P] 

where D is a constant such that D > e 2 . The proposed boundary operator here has the form 


,dt 


+ 6 


du 

dr 


d 1 c > 

dt + 2R + 8 J 


du 

dt 


- 6G{t)u. 


(50) 


We compare this operator with the operator in category 4. When 6 is a constant and G = 0 
one retrieves this operator. The function 8(t) is determined by matching the next asymptotic 
correction in space (i.e., O(jfj) )• This yields 

6 = — . 

- 4R 1 — 2 RG 

The details of these calculations will be found in [3], 

We demonstrate the effect of this proposition through a numerical simulation. For this 
simulation a specific choice of g(t) is considered. It is modeled in such a way that it has 
oscillatory decay. The form is 



u(r 0 ,t) = 1 - 


COS 27T t 

Y+t 2 


t > 0. 


(52) 


Moreover, this problem does not have an explicit analytic solution. Therefore, the numerical 
simulation and comparison is made by doing the calculation in a large domain r 0 < r < R 
with ro = 1 and R = 252. The value of R is chosen in such a way that the waves do not hit 
this artificial boundary so that no reflection in the computational domain is introduced in the 
time period (500 nondimensional units). The specific scheme that is used for computation is an 
explicit second order scheme (both in time and space). The graphs and associated results are 
given in Figure 4. In addition the results that are plotted in the figure represent the magnitude 
of error in time at a fixed point within the computational domain. This point in this case is 
chosen to be r — 2. The results are consistent with the predicted theory. In particular the 
theory predicts that the results of Engquist and Halpern [11] have the largest error for short 
times and is demonstrated in the graph. However, for large times this error tends to decay but 
only at a very short rate which is consistent with the key emphasis of this paper. The second and 
third graphs are due to the computations using the Bayliss and Tuxkel [4] boundary operators 
of orders one and two respectively. Clearly one sees smaller errors at short times and the errors 
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have a tendency of growing larger (they do indeed though not shown here). Finally the last 
curve is the current proposal [3] and has errors less than a percent for all times and it behaves 
well for short, long and in between times. These constructions are a guide to an abbreviated 
proposal that is given for aerodynamic problems, specifically for the Euler equations. 


4 Extensions to Aerodynamics 


Here we extend the ideas presented in the earlier sections to aerodynamic problems governed 
by the Euler equations. A similar analysis can be done for the viscous flows governed by the 
Navier-Stokes equations which is not presented here. It should be noted that to extend the 
analysis of the model problem it is vital to have a linear problem. Clearly the Euler equations 
or the Navier-Stokes equations are not linear. However, the analysis which we seek is for the 
fax field boundary operators and is reasonable to assume linearity of the equations there. The 
form of the Euler equations we use here is as follows: 

^ + div(pu) = 0 (53) 


P 


Du 

P D^ = ~ VP 


De 

Dt 


+ pdivu ~ 0 


(54) 

(55) 


There are a variety of ways to linearize these equations about a state. In this case we consider 
the state at infinity. Since we are dealing with a system of first order equations, additional 
boundary conditions are required according to the hyperbolic theory. The treatments that are 
presented here are for subsonic inflow as well as outflow. Along these lines Gustafsson [18] 
proposed a linearization with the variables q = divu, the divergence of the velocity field, and 
the entropy S = In Noting that u> = curlu the linearized equations look like 


dS dS 

(56) 

a? + “»fe=° 

du> duj 

(57) 

~dt +Uoo ~dx = ° 

S +««,££ + p- ap = ° 

+ u oo2)£ + PooC^oQ = 0 

(58) 
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In the last two equations, elimination of q yields a convective wave equation for the pressure 

(l +u °4) j ,=c ~ A > > - < 59 > 

This is the place where the connection to our earlier analysis holds. The treatment of boundary 
conditions for the pressure is subject to the structure of the solutions of this equation. Before 
we proceed to discuss these equations we point out that the linearized equations for the entropy 
and vorticity suggest that at points of inflow one must specify S and u>. These are dictated by 
the hyperbolic theory as well. 

Now we set up an appropriate initial boundary value problem that will be used to extend 


the analysis in the previous section. The model to be considered here is: 

Tptt + 2 EW*t + Ulrp xx = C'^V'xx + CltPyy (60) 

xj; = ip t = 0 at £ = 0 (61) 

where rj> = p — Poo- Taking the Laplace transform with respect to t yields: 

s 2 i> + 2U 00 s 4> x + UlJ) xx = ClJ) xx + Clrfyy. (62) 

Now we make a dependent variable change from ip to <f> in the form 

$ = e Xx <f>. (63) 

This yields the following calculations: 

^ = Ae^+e^ x (64) 

4> xx = X 2 e Xx 4> + 2Xe Xx <j> x + e Xx <j> xx (65) 

t yy = e Xx K ( 66 ) 

Substituting (64-66) in the transformed equation we obtain 

s a 4> + 2U 00 sX<t> + 2U O0 s<f> x = (Cl - Ul)(<t> xx + 2X<j> x + A V) + C 2 ^. (67) 


We pick the free parameter A (to eliminate the first derivative term) so that 


The remaining terms involving A yield 


s 2 + 2U ooS \ - {Cl - Ul) A 2 = s 2 + 

^ oo ^oo 

a . Ul» 2 * 2 Cl 

ci-ui ci-ui’ 


Thus the equation for (j> becomes 


Cl -VI 


• 2 Cl 
Cl -VI 


<f> = (Cl - Ul)<l> XX + Cl<f>yy 


or 


Let 


'2 «. r /2 ^ — Moc^xx + ^yy* 

oo ^oo 


Vci-vi 

X 

VC-Ml 

y~y- 


Then the equation for <f> becomes 
Radially symmetric solutions satisfy 


^xx H~ fiyy — 5 <f>. 


_ . Kq(st) 
<Pf - S 

K 0 (sr) 


Like in our wave equation calculations for I large we have 




(j + h + f)j 

5 + 5 


<t> 


where 


f — <Ji 


x 2 + y 2 = 


+ y 2 - T 


1 - M to 

r 


' cos 2 6 
l-M r 


— h sin 2 6 




^1 - Mqo sin 2 6. 


Therefore the chain rule yields 


d$ _ s/1 - Ml d$ 

df sjl - Ml sin 2 6 dr ’ 


(69) 

(70) 

(71) 

(72) 

(73) 

(74) 

(75) 

(76) 

(77) 

(78) 

(79) 


15 


Note that 


Therefore 


Noting 


<f> = rpe~ Xz = ip e~ Xrco 80 
d<f> ( dip 


U x s 


— A cos Oip ) e 


— Ar cos0 


A = 




Cl - Ul - Ml 


( 80 ) 

(81) 

(82) 


and letting H{8) — ^ lM °? we see that 

V ’ y/ 1 -Musin' (6) 


Therefore 


H(6) 


dip 


'drP 


(s+ t-H[0) + S) 
ar V TVs 


Sip. 


dr H{6) 


1 (5+ ±H(0) + 6\ Moo . n~j 


dip 

dr 


i + S ) 

'm + h + ««(»)' 


s -f 6 


sip 


where 


J\ - Ml sin 2 6 - Moo cos 6 

m - ' v r-Mi 


(83) 

(84) 

(85) 

( 86 ) 


Upon simplification we propose a boundary operator for the difference of pressure ip = p - poo 
in the form 

(si + { ) *- = - (^1 + n + ss ^) ■ s* - «*<*•*» < 87 > 

where d = y/Cl — Ul. Again one can establish the form of G as 

^ Ri^ +D y m 

Corrections to O ( ) yield in this case 


1 1 


6 = -~= 


(89) 


\Jl - Ml sin 2 6 1 ^RG 

In conclusion, this paper points out potential problems in computing solutions of wave-like 
equations in exterior regions. Remedy of the problems dictates logarithmic correction terms in 
posing radiation conditions. An extension is proposed to the Euler equations in the far field. 
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